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Abstract: Large N QCD is mostly governed by planar diagrams and should show 
linear confinement when these diagrams are suitably summed. The linear confine- 
ment of quarks in a class of these theories using gravity duals that capture the loga- 
rithmic runnings of the coupling constants in the IR and strongly coupled asymptotic 
conformal behavior in the UV was studied in our previous work. We also extended 
the theories to high temperatures and argued the possibilities of meltings and sup- 
pressions of heavy quarkonium states. In this paper we give a formal proof of melting 
using very generic choices of UV completions, and point out some subtleties associ- 
ated with meltings in generic large N theories. Our proof requires only the existence 
of well defined UV behaviors that are devoid of Landau poles and UV divergences of 
the Wilson loops, allowing degrees of freedom to increase monotonously with energy 
scale. We determine the melting temperatures of heavy quarkonium states, which 
could suggest the presence of deconfinement phase transitions in these theories. 
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1. Introduction 

Large N QCD in 3 + 1 dimension is a relatively simpler theory than its finite N 
counterpart, owing in part to the 1/N suppressions of the non-planar diagrams, but 
its solution is still a challenge. For example its not clear how to add up all the 
planar diagrams and argue for linear confinement of fundamental quarks. Although 
in terms of mesonic and glue-ball degrees of freedom the theory may look free, this 
is an oversimplification. Owing to the existence of a Hagedorn temperature and, 
at a more fundamental level, quark degrees of freedom, the system is still complex 
with no simple way of computing, for example, the current-current correlator or the 
master field 

In principle, gauge-gravity duality provides a way to compute - or at least allow 
some analytic control on - some of these quantities. However the problem is that to 
restrict everyting to the supergravity level, where we can have more analytic control, 
the gauge theory should be at strong 'tHooft coupling. This cannot happen for large 
N QCD that is asymptotially free (although there could be a full string theory dual 
description). If one relaxes that condition and looks for strongly coupled asymptotic 
conformal behavior then one can construct a supergravity dual for a class of these 
theories. 

In a previous paper ||, which was a continuation of earlier work ||, a new super- 
gravity background that captures the logarithmic runnings of the coupling constants 
of a particular class of large N QCD in the far IR and the strongly coupled asymp- 
totic conformal behavior in the far UV was constructed. In the intermediate energy 
scales, our dual gravitational background captures the interpolating behavior of the 
beta function. 
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In it was argued that such a geometry would consist of three regions termed 
region 1, 2 and 3 that would capture the IR, the intermediate scale, and the UV 
respectively. The supergravity solution was constructed using fluxes sourced by N 
number of D3 branes, M number of D5 branes and anti-branes respectively while 
taking the back reaction of Nf number of seven branes and Nf — 3 number of anti 
seven branes. The metric in all three regions can be written in the following form |5| 



ds 2 = — = 
Vh 



— g\dt 2 + dx 2 + dy 2 + dz 2 



g 2 g rr dr + g mn dx m dx n (1.1) 



with g,i being the black-hole factors and h being the warp factor depending on all 
the internal coordinates (r, 9i, (pi, ip). To zeroth order in g a Nf and g s M we have our 
usual relations: 

h [0] = 9 [0] = 1 - T f v g& = 1, 9 [ ±dx m dx« = ds 2 Tll 

g l ? = g? = g = l-£, L^g s Na' 2 (1.2) 

where T 11 is the base of a six dimensional conifold and it has the topology of S 3 x S 2 , 
where S 3 and S 2 are three-sphere and two-sphere. Here g s is the string coupling, 
L is the AdS throat radius and is the black hole horizon. The above metric is 
that of a ten dimensional geometry with coordinates (t, x, y, z, r) describing a five 
dimensional non compact space while internal coordinates (9±, 02, 4>i, 02, describe 
a five dimensional compact space. 

At higher orders in g s Nf, g s M, the warp factor, the black hole factors 1 and the 
internal metric get modified because of the back-reactions from the seven-branes, 
three form fluxes and the localized sources we embed. We can write this as: 

h = + fcW, g„ = gf} + g$, g mn = g [ ± + g [ l, 9i = gf ] + gfil.S) 



where the superscripts denote the order of g s Nf and g s M. 

In regions 1 and 2, the warp factor h has logrithmic dependence while in region 
3, it can be expanded exclusively as a power series in 1/r = u. In terms u coordinate, 
we have the following Taylor series expansion for the warp factor h: 

-L = A n u n - 2 (1.4) 

where A n are Taylor coefficients and in general can be functions of internal coordi- 
nates. Region 3 was taken to be large enough so that the Nambu-Goto string would 

1 This also means that the black hole horizon is no longer the surface r — and the equation 
for the surface may be more involved. 
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lie completely inside it. The computation of the potential at zero temperature re- 
vealed that at short distances the potential should be dominated by the Coulomb 
term 2 !!: 



Vqq = ~ — — + 0{d) (1.5) 



0.236 
~d 

where the numerical value and the sign of the first term are fixed naturally by the 
dual background. It is interesting to note that this value is of the same order of 
magnitude as the one derived from QCD lattice simulations ||. For large distances 
one expects the potential to be dominated by the linear term ||: 

A T n 



V QQ = ( ) d (1.6) 

" A max 



where x max is the maximum value for the depth of the U shaped string denoted by 
M max and it is a solution of the constraint equation (see section 2 and also || for 
details): 

i(m + l)An+ 3 < a + x 3 = 1 (1.7) 

The constraint equation is obtained by demanding that the distance between the 
quarks and the potential energy of the quark pair is real. For AdS space, interquark 
distance d and potential Vqq is always real for any value of w max , and hence x max do 
not exist, as there in no constraining equation. 

At high temperatures and density, medium effects should screen the interaction 
between the heavy quark the anti-quark. The resulting effective potential between 
the quark anti-quark pairs separated by a distance d at temperature T can then be 
expressed succinctly in terms of the free energy F(d,T), which generically takes the 
following form: 

F(d,T) = adf s (d,T)-^f c (d,T) (1.8) 

where a is the string tension, a is the gauge coupling and f c and f s are the screening 
functions 3 (see for example || and references therein). For the quark and the anti- 
quark pair kept at respective positions of +~ and — |, we expect the Wilson line 
W to be related to the free energy through 



exp 



F(d,T) 



(w* (+I))(w(-D) 



2 In deriving the Coulomb potential, we have set AdS throat radius L = 1 and string tension 
Tq = 1 for convenience, and this is only a redefinition of units. By restoring units, one obtains 
Vqq = - a236 / 57W + 0(d), as expected [@. 

3 We expect the screening functions f s , f c to equal identity when the temperature goes to zero, 
yielding the zero temperature Cornell potential. 
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In terms of Wilson loop, the free energy Eq.(L8) is now related to the renormalised 
Nambu-Goto action for the string on a background with a black- hole 4 . One may 
also note that the final theory is not three-dimensional, but four-dimensional and 
compactified on a circle in Euclideanised version. 

Using Eq. (|1.9|) and the identification of the Wilson loop to the Nambu-Goto 
action, the free energy (or equivalently the potential) between the heavy quark and 
the anti-quark at non-zero temperatures can be deduced. For large distances the 
potential is [|j 



v ««^i 1 - y -k{^b) d (L10) 

where y max is the maximum value for the depth M max of the U shaped string in the 
presence of a black hole. It is given by solving the constraint equation: 



\{m + l)A m+3 y^ + I {J (k - I) ' 3 ' 

J' h—n \ ' 



Ay l max ( \ + 2 j - i 



which arises by demanding that inter quark separation and free energy is always 
real ||. Note that Uh = 1/r/j and for zero temperature, Uh = oo. Thus at zero 
temperature, Uh — > oo and ( |1.11| ) reduces to (|1.7|). 



At this point, it is important to define in a more precise manner what exactly is 
meant by melting. In this work, "melting" is meant to quantify the disappearance 
of the linear portion of the quark- ant iquark potential, at a given length scale. There 
is therefore a qualitative difference between this behavior and the actual dissolution 
of quarkonium bound states: The robustness of the QQ pairs will actually depend 
on where their energy sits in the potential profile. Later studies could also involve 
a more precise characterization through a study of the temperature dependence of 
quarkonium spectral densities, for example. 

Going back to Eq.( |1.10|) , one might think that the melting temperature u~ x is 



given by the condition that y max = Uh is a solution of Eq. ( |1.11[ ) so that the coefficient 
of d in Eq. flL.lO ) vanishes. Actually, this turned out not to be the case. As shown 



below, Eq. fll.llD can never have a solution at y m3X = u^. Hence, the deconfinement in 
our case does not mean a complete absence of the linear potential. Rather, it means 
that the linear potential ceases to have an infinite range and as temperature increases, 
the range of the linear potential quickly becomes short provided that Eq. ( p..ll| ) allows 



4 Thcre is a large body of literature on the subject where quark anti-quark potential has been 
computed using various different approaches like 0. Although the potential Eq.(LfO) that we get 



matches well with other results, we have an additional constraint: Eq.( f .1 1 ), from the background 
RG flow. This additional constraint, which cannot be seen from an AdS or an AdS with an IR 
cut-off dual, will have non-trivial consequences for melting that we will discuss in section 2. 
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a real positive solution lying in region 3. It is not clear this would always be the 
case, so Eq . (|1 . 1 1| ) requires a more detailed investigation. In the following section i.e 
sec. 2, we will carefully analyse Eq. ( |1.11| ) for generic choices of warp factors or more 
appropriately, generic UV completions that have no Landau poles or UV divergences 
of the Wilson loops, and give a proof of quarkonium melting for this class of theories. 
Section 3 contains a detailed numerical analysis that will allow us to find the melting 
points, and could also suggest the presence of a deconfinement phase transition in 
these theories. We conclude with a short discussion. 



2. Proof of the existence of a melting point 

Our cascading picture of renormalization group flow demands that in the region 3, 
the effective number of colors grows as r grows. The number of colors at any scale 
u — 1/r in the region 3 is given by N e g(u) = N(l + aiu 1 ). For the analysis given 
here, it is simpler to define and use 



T(u) = ^-= lK 1 ( 2 1) 

instead of N e s(u). The coefficients A n appearing in Eqs. (|1.4j), ( |1.7| ) and ( |1.11| ) are 
related to F{u) by J-'iu) = A n u n ; and h is the warp factor. In terms of Fiu), the 
condition that N eS (u) is a decreasing function of u = 1/r becomes 

F'(u) > (2.2) 

Combining Eqs. (|2. 1|) and (|2.2| ) yields the following condition 

T(u) > ^ (2.3) 
From now on, and as mentioned earlier, the value of L is set to 1 so that J-'iu) > 1. 
2.1 Zero temperature 

Let w m ax be the maximum value of u for the string between the quark and the anti- 
quark. Then the relationship between w max and the distance between the quark and 
the anti-quark is given by 



"(Wmax) = 2M max y r (M max ) / dV 



(^(Wmax«)) 2 



.,-1/2 



JF{u v 



(2.4) 



Two conditions must be met before asserting that this expression represents the phys- 
ical distance between a quark and an anti-quark in vacuum. One obvious condition 
is that the integral is real. This is guaranteed if for all < v < 1: 

W(v\u max ) = v 2 ( ^ Um3X \ ) < 1 (2.5) 
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Another condition is that the potential between the quark and the anti-quark must 
be long ranged. That is, d(u max ) must range from to oo as u m3X varies from to 
some finite value, say M max = x max . Since J-{u) > 1, the only way to satisfy these 
conditions is via the (sufficiently fast) vanishing of the square-root in Eq. ( |2.4|) as 

V y 1 at Umax 3?max- 

For most w max , 1 — W \v\u max ) 2 vanishes only linearly as v approaches 1. In this 
case, d(u max ) is finite since the singularity in the integrand behaves like 1/a/1 — v 
and hence integrable. To make <i(w max ) diverge at w max = x max , 1 — W(v |x max ) 2 must 
vanish quadratically as v approaches 1 so that the integrand is sufficiently singular, 
1/a/1 — W(v\x max ) 2 ~ 1/|1 — v\. Therefore, the function W(v\x max ) must have a 
maximum at v — 1. 

To determine the value of x max , consider 

T rW/ I \ o I -^""(-^max) \ / -, / \ J~ {X max V) \ /~ n \ 

W'(v\x max ) = 2v — y - J - 1 - (x max v)—± !r (2.6) 

K-Fixm^v) J V 2J r (x max V) J 

For this to vanish at v — 1, a; max must be the smallest positive solution of 

xF\x) - 2JF{x) = (2.7) 

With the definition ^(u) = A n u n , one can easily show that this is equivalent to the 
condition ( |1.7|) which was originally derived in ||. The allowed range of M max is then 

^max — 3<max (2.8) 

and within this range, d(u max ) varies from to oo. How it varies will depend on the 
values of Q m as well as T{u). 

2.2 Finite temperature 

At finite temperature, the relation between w max and the distance between the quark 
and the anti-quark is obtained by replacing T{u) with ■sjx — u 4 /uf l J r (u) in Eq.([2.4|): 



d T (u mBX ) = 2u max Wl - Umax/^fo-^max) / dv 



1 „,2 



X _ V A (1 - U tJ U h) ( ^(«max) X ' 



i-vW m jui)(F( Umax v)y 

-1/2 

(i-^ 4 <axK) j (2,9) 

The explicit factor of w max makes dT(u max ) vanish at u max = as in the T = case. 
As M max approaches Uh, the integral near v = 1 behaves like 

rfr W ~ fdv ^~<J< (2.10) 
Jo a/(1 - v)(l - vu max /u h ) 

which indicates that dT{u max ) goes to as u max approaches u^. Hence, at both M max = 
and u max = Uh, dT(u max ) vanishes. Since dT(u max ) is positive in general, there has 
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to be a maximum between u max = and w max = Uh- Whether the maximum value of 
^T(wmax) is infinite as in the T = case depends on the temperature (equivalently, 
u h l ) as we now show. 

The fact that the physical distance needs to be real yields the following condition. 
For all < v < 1, 



HUm^v) J V 1 - ui lax V 4 /u h 



w T (v\ Umax ) = max : j- — < i (2.ii) 



Taking the derivative gives 

(1 - Wmax/?4) 1/2 ^(Mmax) 



Wt( U l M max) 



(l-<ax^ 4 K) 3/2 ^max^) 

X [-(Mmax^)(l - (u max V /Uh) 4 )^' (U max v) + 2J r (u max v)] 



(2.12) 



Similarly to the T = case, d,T(u max ) can have an infinite range if the derivative 
vanishes at v — 1 for a certain value of M max , say M max = 2/max- This value of y max is 
determined by the smallest positive solution of the following equation 

yT\y) - 2F(y) = {y/u h ) 4 yF{y) (2.13) 

which then forces Wf(l|y max ) to vanish. Note that the left hand side is the same 
as the zero temperature condition, Eq. fl2.7|) . The right hand side is the temperature 
(uh) dependent part. Using the facts that: 

j—l oo 

H {k - 1/2) = -(2j - 3)\\/2\ x3 & ~ : ^ !! = 1 " v / T^, (2-14) 

k=0 7=1 



it can be readily shown that Eq.( p.l3|) is equivalent to Eq . ( |1 . 1 1|) as long as M max < uu 



It is also clear that y = u h cannot be a solution of Eq. (|2.13| ) because at y — Uh, the 
equation reduces to J-(uh) = which is inconsistent with the fact that J-"(y) > 1. 

Recall that JF{y) > 1 and J^'iy) > in the region 3 and we assume that the 
equation yF'(y) — 2J-"(?/) = has a real positive solution x max . Hence as y increases 
from towards x max , the left hand side of Eq.( |2.13| ) increases from —2 while the right 
hand side increases from 0. The left hand side reaches when y = x max which is the 
point where the distance rf(u max ) at T = becomes infinite. At this point the right 
hand side of Eq. (|2.13| ) is positive and has the value (x max / , u ?i ) 4 x max J : ''(x max ). Hence 
the solution of Eq. ([2.13|) , if it exists, must be larger than x max . 



Consider first low enough temperatures so that Uh ^> x max . For these low tem- 
peratures, Eq. ( |2.13| ) will have a solution, as the right hand side will be still small 



around y = x max . This then implies that the linear potential at low temperature will 
have an infinite range if the zero temperature potential has an infinite range. 
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To show that the infinite range potential cannot be maintained at all tempera- 
tures, let Uh = x max . When the left hand side vanishes at y = x max , the right hand 
side is x max J-"'(x max ) = 2J r (a; max ) which is positive and finite. For y > x max , the left 
hand side (yF'iy) — 2J r {y)) may become postive, but it is always smaller than yF'ii/) 
since J^iy) is always positive. But for the same y, the right hand side ( y (y/uh) 4 yJ-''(y)) 
is always positive and necessarily larger than yF'^y) since (y/uh) > 1. Hence, 
Eq.( |2.13 ) cannot have a real and positive solution when Uh = x max . Therefore be- 



tween Uh = oo and Uh = i ma x, there must be a point when Eq. (|2.13|) cease to have a 
positive solution. 

When Eq. (|2.13|) has no solution, then the expression for (^(wmax), ( |2.9|) will not 
diverge for any u max within (0,Uh)- Furthermore, since the expression vanishes at 
both ends, there must be a maximum cl T (« max ) at a non-zero u max . When the distance 
between the quark and the anti-quark is greater than this maximum distance, there 
can no longer be a string connecting the quark and the anti-quark. 

3. Numerical analyses of melting temperatures 

After discussing the most general choice for warp factors that give rise to y max and 
consequently linear potential, specific examples of geometries that may arise as so- 
lutions to Einstein's equation will be considered, starting with the following ansatz 
for the metric: 

ds 2 = — ^=dt 2 + -^=(dx 2 + dy 2 + dz 2 ) + ( -^du 2 + ds 2 M 
Vh Vh u 2 \gu 2 

- 9 -.dt 2 + -^=(dx 2 + dy 2 + dz 2 ) + ^-g mn dx m dx n (3.1) 



where h = h(u, 8i, fa, ip),H = H(u, 9i, (pi, ip), and g = 1 — w 4 /w^; is the compact 
five dimensional manifold parametrised by coordinates (9i, (pi, ip) and can be thought 
of as a perturbation over T 1,1 . Here u = is the boundary and u = Uh is the 
horizon. As discussed in earlier work ||, the above metric arises in region 3 of |§ 
when one considers the running of axio-dilaton r, D7 brane local action and fluxes 
due to anti five-branes on a geometry that deviates from the IR OKS-BH geometry 
from the backreactions of the above sources. The three-form fluxes sourced by (p, q) 
anti-branes are proportinal to r~ l f(r) for some positive i (see for details about 
/(r)), where the function f(r) — > 1 as r — > oo and f(r) — > as r — > 0. With the 
coordinate u = 1/r, there is another function: k{u) = exp(— u A ), A > 0, that also 
has somewhat similar behaviour as f(u) and may allow us to have a better analytic 
control on the background. With such a choice of k(u), the total three from flux is 
proportional to u A M{u) with 

M(u) = Af [1 - k(u)) = M [l - exp(-u A )] (3.2) 
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where M is the number of bi-fundamental flavors. Thus three-form fluxes are de- 
caying fast as Mu A [l — exp (— u A )~\ and, as shown in ||, the seven-branes could be 
arranged such that the axio-dilaton r behaves typically as r ~ u B . This means that 
from the behavior of the internal Riemann tensor one may conclude that the internal 
metric g mn behaves as g mn ~ tt c exp(c w c ) where A,C,A and C are all positive and 
c Q could be positive or negative depending on the precise background informations. 

From the above discussions it should be clear that taking the three-forms and 
world-volume gauge fluxes to be exponentially decaying in the IR (but axio-dilaton 
to be suppressed only as u B ) should solve all the equations of motion, giving the 
following behavior for the warp factor h and the internal metric H in (|3.1|) 5 : 

h = LVexp(-cra 5 ), H = exp(ftj) (3.3) 

where we are taking a, a, ft, ft to be all positives with a, ft to be functions of internal 
coordinates (6*j, <pi, ip) and L 4 = g s Na' 2 to be the asymptotic AdS throat radius 6 . 

Motivated by the above arguments, we will consider Nambu-Goto action of the 
string in the geometry with (a, (3) = (3,3) and (5, ft) = (4,4) at temperatures 
and respectively in Eq.( [3.3|) . As in [[| we consider mappings X^(a, r), which 
are points in the internal space, to lie on the slice: 

6 1 = 9 2 = vr, 0, = 0, ^ = (3.4) 

so that on this slice a, ft are fixed and we set it to (a, ft) = (0.1, 0.05) for both choices 
(5, ft). (Such a choice of slice will also help us to ignore the three-form contributions 
to the Wilson loop.) With these fixed choices for the warp factors, we plot the 
interquark separation d as a function of u max i n figures 1 and 2 for various values of 
T = l/u h . 

Note that for both choices of warp factors, for low enough temperatures, there 
exist w ma x = Z/max where d —¥ oo. As the temperature is increased, y max increases 
modestly. On the other hand from figure ||, one sees that when T > ~ 0.28 
there exists a d m8/X which is finite. This means for interquark distance d > d mgx , there 
is no string configuration with boundary condition x(0) = ±d/2 implying that the 
string attaching the quarks breaks and we have two free partons for d > d max . Thus 
we can interpret d max to be a "screening length" . From figure 2 we observe similar 

(2) 

behaviour but now d max exists for T > T c ~ 0.399. 

In figure 3, d max as a function of T is plotted. We note that for a small change 
in the temperature near (or near equivalently) there is a sharp decrease in 
screening length <i max , but for T » Tc \ % = 1,2, the screening length does not 



5 See also the interesting works of || where exponential warp factors have been chosen. 

6 Note that /? in (3.3) could be considered negative so that H would be decaying to zero in the IR. 



However since region 3 doesn't extend to the IR we don't have to worry about the far IR behavior 
of Eq.O). 
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Figure 1: Interquark distance as a function of tt max f° r various temperatures and warp 
factor with (a, a, ft, ft) = (0.1,3,0.05,3) in the warp factor equation. 




Figure 2: Quark-antiquark distance as a function of u mSuX for various temperatures and 
warp factor with (a, a, ft, ft) = (0.1,4,0.05,4) in the warp factor equation. 



change much. In fact c/ max behaves as C + exp(— jT) (where C and 7 are constants) 
which in turn could be an indicative of a phase transition near T c (i) fori = 1,2 i.e 
the two choices of warp factor. 

Finally we plot the potential Vqq as a function of d in figures 4 and 5 for the 
two choices of warp factor. For T < in figure 4 and T < in figure 5, we 
have potential energies linearly increasing with an arbritrarily large increment of the 
interquark separations. Thus we have linear confinement of quarks for large distances 
and small enough temperatures. For T > Tc % \ i = 1 or 2, there exists a d max and 
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Figure 4: Heavy quark potential Vqq as a function of quark separation d with cubic 
warp factor, or equivalently, (a,a,(3,/3) = (0.1,3,0.05,3) in the warp factor equation for 
various temperatures. As mentioned in the text, one shouldn't consider the free energy 
(or equivalently the potential) to stop abruptly in the plot. After the string connecting 
the quarks breaks, the curves should be extrapolated by the contributions to the potential 

(i) 

energy from the two free strings and their world- volume fluctuations for all T > T c ■ . 



for all distances d > <i max there are no Nambu-Goto actions, Sng, for the string 
attaching both the quarks. This means that we have free quarks and Vqq is constant 
for d > d max . Of course looking at figures 4 and 5 one shouldn't conclude that the free 
energy stops abruptly. What happens for those two cases is that the string joining 
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Figure 5: Heavy quark potential Vqq as a function of quark separation d with quartic 
warp factor, or equivalently, (a, a, /3,/3) = (0.1,4,0.05,4) in the warp factor equation for 
various temperatures. Here again, beyond T > T c the curves should be extrapolated by 
the contributions to the potential energy from the two free strings and their world- volume 
fluctuations. 



the quarks breaks, and then the free energy is given by the sum of the energies of 
the two strings (from the tips of the seven-branes to the black-hole horizon) and the 
total energies of the small fluctuations on the world-volume of the strings. The latter 
contributions are non-trivial to compute and we will not address these in any more 
detail here, but energy conservations should tell us how to extrapolate the curves in 
figures 4 and 5, beyond the points where the string breaks, for all T > T c . Of course 
after sufficiently long time the two strings would dissipate their energies associated 
with their world-volume fluctuations and settle down to their lowest energy states. 

Observe that for a wide range of temperatures < T < T c , the potential 
and thus the free energy hardly changes. But near a narrow range of temperatures 
T c — e < T < Tc + e (where e ~ 0.05), free energy changes significantly. For figure 
5 the change is more abrupt than figure 4. This means as we go for bigger values of 
5, the change in free energy is sharper. 

In figure 6, we plot the slope of the linear potential as a function of T. Again 
for a wide range of temperatures, there are no significant changes in the slope but 
near Tc % \ the change is more dramatic: the slope decreases sharply, indicating again 
the possibility of a phase transition near Tc \ As we noticed before, here too bigger 
exponent 5 gives a sharper decline in the slope hinting that when a » 1, the 
transition would be more manifest. 

To conclude, the above numerical analyses suggest the presence of a deconfine- 
ment transition, where for a narrow range of temperatures 0.28 < T c < 0.39 the free 
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Figure 6: Slope of linear potential as a function of T for both cubic and quartic warp 
factors. Note that in the figure the slopes have been computed as with a range of d 
from 1.6 to 1.7 in appropriate units. 



energy of QQ pair shows a sharp decline. Interestingly, changing the powers of u 
in the exponential changes the range of T c only by a small amount. So effectively 
T c lies in the range 0.2 < T c < 0.4. Putting back units, and defining the boundary 
temperature 7 T as T = 9 , our analyses reveal: 

4ny/h(u h ) 

0.91 1.06 , , 

~ix - Tc ~ (3 ' 5) 

which is the range of the melting temperatures in these class of theories for heavy 
quarkonium states. Since the temperatures at both ends do not differ very much, 
this tells us that the melting temperature is inversely related to the asymptotic AdS 
radius in large N thermal QCD. 



4. Discussions and conclusions 

In this work we attempted a small step towards the understanding of the melting 
temperatures of heavy quarkonium states in a class of large iV thermal QCD that 
have well defined UV behaviors without Landau Poles and UV divergences of the 
Wilson loops. Our analyses reveal two interconnecting results: any large iV gauge 
theory with Nf fundamental quarks (Nf N) that is away from conformality with 
atmost one UV fixed point 8 should always have a mass gap and, consequently at 

7 See sec. (3.1) of || for details. 

8 This also means a good UV behavior with degrees of freedom increasing monotonously with 
energy scales. 
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a certain temperature T c , the heaviest quarkonium states in such a theory should 
show melting. Such a conclusion seems to be true for generic non-supersymmetric 
theories, and it would be interesting to extend this to the supersymmetric cases. 

The melting temperatures 9 that we get for various UV completions seem not 
to differ too much from each other. This may mean that there is some underlying 
universal behavior of heavy quarkonium states in large N theories with good UV 
behaviors at strong 'tHooft couplings. An interesting related study would be that of 
the behavior at weak 'tHooft couplings. There, there are also Hagedorn states that 
would come in because of the unsuppressed string modes. Details on this will be 
presented elsewhere. 
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